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Abstract. The gap equation is a cornerstone in understanding dynamical chiral symmetry breaking and 
may also provide clues to confinement. A symmetry-preserving truncation of its kernel enables proofs of 
important results and the development of an efficacious phenomenology. We describe a model of the kernel 
that yields: a momentum-dependent dressed-quark propagator in fair agreement with quenched lattice- 
QCD results; and chiral limit values: f° = 68MeV and (qq) = — (190 MeV) 3 . It is compared with models 
inferred from studies of the gauge sector. 

PACS. 12. 38. Aw General properties of QCD - 11.30.Rd Chiral symmetries 



1 Introduction 

We will address these topics from the perspective of QCD's 
Dyson-Schwinger equations (DSEs) Q. The DSEs are a 
keystone in proving renormalisability and provide a gen- 
erating tool for perturbation theory. The latter point is 
very important in applications to low-energy phenomena 
because it means that the mo del- dependence which nec- 
essarily appears in continuum studies of nonperturbative 
QCD is restricted to the infrared; i.e, to momentum-scales 
~ 1 GeV 2 . This feature has successfully been exploited in 
applications to the spectra H,||, and strong |Q] and elec- 
troweak [|| interactions of mesons. It also mitigates, to 
a useful extent, some of the problems with the approach; 
e.g., it helps in developing reliable truncations for the cou- 
pled system of DSEs. 

Recent successes are founded on an accurate under- 
standing of dynamical chiral symmetry breaking (DCSB), 
and its role in resolving the dichotomy of the pion [|| (as 
both a Goldstone mode and a massless bound state of 
massive constituents) and its relation to the long-range 
behaviour of the effective coupling between quarks . 



2 Dynamical chiral symmetry breaking 

This is a purely nonperturbative phenomenon that can be 
studied via QCD's gap equation: 



wherein: m is the current-quark bare mass; g is the cou- 
pling constant; D^ v {p—q) is the dressed-gluon propagator; 
r®(q,p) is the dressed-quark-gluon vertex; and the solu- 
tion is the dressed-quark propagator 



S(p) 



Z{P 2 ) 



ij ■ pA(p 2 ) + B(p 2 ) ij-p+M(p 2 ) 



(2) 



(Equation (|l]) is the unrenormalised equation: renormali- 
sation will only be mentioned as necessary.) 

As noted in the Introduction, one can use the gap equa- 
tion to evaluate the fermion self-energy perturbatively and 
thereby obtain 
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S(p) 1 — (*7 ' P + m ) 
d 4 q 



(27T) 



g 2 D^{p- q )—^S(q)r-(q,p), (1) 



wherein the ellipsis denotes two- and higher-loop contri- 
butions. It is apparent that each term is proportional to 
the current-quark mass and vanishes as m — > 0. Hence, 
in perturbation theory, if there is no quark mass to be- 
gin with, interactions do not generate one; i.e., DCSB is 
impossible in perturbation theory. 

In the chiral limit, namely, m = 0, Eq. (Q) always 
admits the trivial solution B(p 2 ) = 0; i.e., the solution 
accessible perturbatively. However, as is characteristic of 
gap equations, when the integrated strength of the ker- 
nel exceeds some critical value, a B{p 2 ) ^ solution 
is generated. This nonperturbative, dynamical generation 
of a quark mass from nothing is DCSB. The integrated 
strength of the kernel can be characterised by an interac- 
tion tension, a A , for which the critical value is M a A ~ 
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2.5 GeV/fm: this strength provides for DCSB but onl 
just. An acceptable description of hadrons requires 
g a ~ 25GeV/fm, which is an order of magnitude larger. 



2.1 Symmetry preserving truncation 

It is apparent that the gap equation's kernel is formed 
from a product of the dressed-gluon propagator and dres- 
sed-quark-gluon vertex. It may be calculated in pertur- 
bation theory but that is inadequate for the study of in- 
trinsically nonperturbative phenomena. Consequently, to 
make sensible statements about DCSB, one must employ 
an alternative systematic and chiral symmetry preserving 
truncation scheme. 

One must also focus on more than just the gap equa- 
tion's kernel. Chiral symmetry is expressed via the axial- 
vector Ward-Takahashi identity: 



P M r 5fl (k; P) = S-\k+) ^75 + i 7 5 



(4) 



k± = k ± P/2, wherein /5 M (fc;P) is the dressed axial- 
vector vertex. This three-point function satisfies an inho- 
mogeneous Bethe-Salpeter equation: 



d 4 q 

W) 



4 [S(q + )r^(q; P)S(q^)] sr K™(q, k;P), (5) 



in which K(q, k\ P) is the fully-amputated quark-antiquark 
scattering kernel. The Ward-Takahashi identity, Eq. (|J), 
means that the kernel in the gap equation and that in 
the Bethe-Salpeter equation are intimately related. There- 
fore a qualitatively reliable understanding of chiral sym- 
metry and its dynamical breaking can only be obtained 
using a truncation scheme that preserves this relation, and 
hence guarantees Eq. (||) without a fine-tuning of model- 
dependent parameters. 

One such scheme exists ||. Its leading-order term is 
the renormalisation-group-improvcd rainbow-ladder trun- 
cation and the general procedure provides a means to 
identify, a priori, those channels in which that truncation 
is likely to be accurate. This scheme underlies the suc- 
cessful application of a rainbow-ladder model to flavour- 
nonsinglet pseudoscalar mesons and vector mesons 

The systematic nature of the scheme has also made 
possible a proof of Goldstone's theorem in QCD Q and 
the derivation of a mass formula that unifies the light- and 
heavy-quark sectors . 

In this scheme, as in perturbation theory, it is impossi- 
ble, in general, to obtain complete closed-form expressions 
for the kernels of the gap and Bethe-Salpeter equations. 
However, for the planar dressed-quark-gluon vertex de- 
picted in Fig. [j], closed forms can be obtained and a num- 
ber of significant features illustrated jll] when one uses 
the following model for the dressed-gluon line fl2|| 




Fig. 1. Planar dressed-quark-gluon vertex obtained by 
neglecting contributions associated with explicit gluon self- 
interactions. Solid circles indicate fully dressed propagators. 
The quark-gluon vertices are not dressed. 



where Q, measured in GeV, sets the model's mass-scale. 
This form is a precursor to the class of renormalisation- 
group-improved models. It has many positive features in 
common with that class and, furthermore, its particular 
momentum-dependence works to advantage in reducing 
integral equations to character-preserving algebraic equa- 
tions. Naturally, there is a drawback: the simple momen- 
tum dependence also leads to some model-dependent arte- 
facts, but they are easily identified and hence not cause 
for concern. 

It is a general result that with any vertex whose 
diagrammatic content is known explicitly; e.g., Fig. [j], it is 
always possible to construct a unique Bethe-Salpeter ker- 
nel which ensures the Ward-Takahashi identities are au- 
tomatically fulfilled: that kernel is necessarily non-planar. 
This becomes transparent with the model in Eq. (^) , using 
which the gap equation obtained with the vertex depicted 
in Fig. |l| reduces to an algebraic equation, irrespective of 
the number of dressed-gluon rungs that are retained, and 
the same is true of the Bethe-Salpeter equations in every 
channel: pseudoscalar, vector, etc. 

Results for the ir and p are illustrated in Table [|. It is 
evident that, irrespective of the order of the truncation; 
i.e., the value of n, the number of dressed gluon rungs in 
the quark-gluon vertex, the pion is massless in the chiral 
limit. (NB. This pion is composed of heavy dressed-quarks, 
as is evident in the calculated scale of the dynamically 
generated dressed-quark mass function: M(0) « 0.5 GeV.) 
The masslessness of the tt is a model-independent conse- 
quence of the consistency between the Bethe-Salpeter ker- 
nel and the kernel in the gap equation. Furthermore, the 
bulk of the p-n mass splitting is present in the chiral limit 
and with the simplest (n — 0; i.e., rainbow-ladder) ker- 
nel, which makes plain that this mass difference is driven 
by the DCSB mechanism: it is not the result of a finely 



Table 1. n and p meson masses obtained with Q = 0.48 GeV. 
(Dimensioned quantities in GeV.) n is the number of dressed- 
gluon rungs retained in the planar vertex, see Fig. [j], and hence 
the order of the vertex-consistent Bethe-Salpeter kernel. 
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adjusted hyperfine interaction. Finally, the quantitative 
effect of improving on the rainbow-ladder truncation; i.e., 
including more dressed-gluon rungs in the gap equation's 
kernel and consistently improving the kernel in the Bethe- 
Salpeter kernel, is a 10% correction to the vector meson 
mass. Simply including the first correction (n = 1; i.e., 
retaining the first two diagrams in Fig. |l|) yields a vector 
meson mass that differs from the fully resummed result by 
~ 1%. The rainbow-ladder truncation is clearly accurate 
in these channels. 



2.2 Rainbow-ladder truncation 

The renormalisation-group- improved rainbow- ladder trun- 
cation is based on the fact that in QCD, on the domain 



for which Q := (fe — q)* 



q is large and spacelike, 



M(q,k;P) :=g 2 (fi 2 )D^(k-q) 

x [r«(k + ,q + )S(q+)] ® [S(g_)r»( g _,*_)] 



A^a(Q z )Dll%k~q) 



X" 



S 



free 



(<?-) 



(7) 



free 



where a(Q 2 ) is the strong running coupling and, e.g., S 
is the free quark propagator. The dressed-ladder trunca- 
tion supposes that Eq. (0) is valid for all momenta and is 
thus an assumption about the long-range (Q 2 ~ 1 GeV 2 ) 
behaviour of the interaction. Requiring that Ward-Taka- 
hashi identities be automatically fulfilled by any trunca- 
tion entails (9[|ll| that the kernel in the renormalised gap 
equation is expressed through 



A" 



G{Q*)D*?{Q)—^S{q)- lv , 



(8) 



where we have introduced G(Q 2 ), an effective interaction, 
to emphasise that, while Q(Q 2 ) — 47ra(Q 2 ) is well-deter- 
mined for Q 2 > 1 GeV 2 , the behaviour at Q 2 < I GeV 2 ; 
i.e., at infrared length-scales (~ 0.2 fm), of the kernels in 
the gap and Bethe-Salpeter equations is unknown. 

Contemporary DSE studies employ a model for the 
infrared behaviour. The most extensively applied is H: 



G(Q 2 ) 
Q 2 

+ 8tt 



hi 



Q 2 M 2 



QCD 



HQ 2 ) 



(9) 



where: f{Q 2 ) = [l-exp(-Q 2 /[4m: 

„ J2 



A. 



(4) 



1 })}/Q 2 ,m t = 0.5 GeV; 
N f = 4; and tIqcd = 
i gives a{m 2 z ) = 0.126. 



+0.009 
-0.013' 



means 



- 1; 7m = 12/(33 - 2N f ) 
vs 0.234 GeV Jl3) (NB. Eq. (| 

Comparison with a modern value ]14| |: 0.113 
that a smaller Aqct> is acceptable in the model, if one 
wants to avoid overestimating the coupling in the ultravi- 
olet, but not a larger value.) The true parameters in Eq. 
(^) are D and w, which together determine the interaction 



tension in the model. However, they are not independent: 
in fitting, a change in one is compensated by altering the 
other; e.g., on the domain u> G [0.3, 0.5] GeV, the fitted 
observables are approximately constant along the trajec- 
tory us D — (0.72 GeV) 3 . This correlation: a reduction in 
D compensating for an increase in to, acts to keep a fixed 
value of the interaction tension; i.e., the interaction's in- 
tegrated strength on the infrared domain. With 



D = (0.96 GeV) 2 , uj 
and current-quark masses 

ml GeV = 5.5 MeV, ml GcV 



0.4 GeV 



(10) 



f 25 MeV, (11) 



one obtains the excellent description of meson observables 
described in the Introduction, which can be illustrated via 
the calculated values of the leptonic decay constant and 
vacuum quark condensate in Table ^. 

2.3 Comparison with lattice-QCD simulations 

The scalar functions characterising the renormalised dres- 
sed-quark propagator: the wave function renormalisation, 
Z(p 2 ), and mass function, M (p 2 ), obtained by solving the 
renormalised gap equation using Eq. (||), are depicted in 
Figs. [| |. The infrared suppression of Z(p 2 ) and enhance- 
ment of M ( p 2 ) are long-standing predictions of DSE stud- 
ies, as Rcf. Mil makes plain and could be anticipated from 
Ref. This prediction has recently been confirmed in 
numerical simulations of quenched lattice-QCD, as is evi- 
dent in the figures. 

It is not yet possible to reliably determine the chiral 
limit behaviour of lattice Schwinger functions. Hence one 
does not have a lattice estimate for the quantities in Table 
||. To obtain such an estimate, we used the DSE model 
described previously, varying (D,u>) in order to reproduce 
the lattice data. Our current best result, obtained with 



D = (0.74 GeV) 2 , lu = 0.3 GeV , 



(12) 



at a current-quark mass of 0.6m] GcV , Eq. (pi]), is also 
depicted in Figs. ||, [|, and yields = 0.094 GeV, — 
0.48 GeV. (NB. These are values for a pion-like bound 
state constructed from constituents whose current-quark 
mass is to « 14 m„. Our DSE-calculated meson-mass ver- 
sus current-quark-mass trajectory is consistent with re- 
sults of recent quenched lattice-QCD simulations obtained 
on m G [1,2]to s [ffoj.) The parameters in Eq. ( p"2| ) give 
chiral limit results: 



/° = 0.068 GeV, 



-(9<?}lGcV 



(0.19GeV) E 



(13) 



Table 2. Calculated pion decay constant and vacuum quark 
condensate. Experimentally, f w — 0.092, and a best plie- 
nomenological value of the condensate is : — (gg)iGeV — 



(0.236 ± 0.008) d 



(All quantities are listed in GeV. 



Model 



Eqs. (|L (M) 
Eqs. (14), (15) 
Eqs. (~), (16) 



{-mi Gc v) iri 



0.092 
0.053 
0.099 



0.24 
0.f5 
0.25 
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DSE soln, Eqs. (9) and (12) 
lattice Z(p) 
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Fig. 2. Wave function renormalisation. Solid line: solution 
of the renormalised gap equation using Eqs. (^|), ([[o|); data: 
lattice simulations 1161 , m — 0.036/a ~ 60MeV; dashed- line: 
gap equation solution using Eqs. (^), (p^). (In the DSE studies, 
the renormalisation point £ = 19 GeV and the current-quark 



mass is 0.6 m 



1 GeV 



[Eq. (hjM) 



0.5 - s 



0.4 - \ 



I 1 I 1 I 1 I r 
DSE soln, Eqs. (9) and (10) 
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lattice M(p) 



\ 

I \ 



0.3 



0.2- 



0.1- 



0.5 



1.5 2 2.5 

p (GeV) 



3.5 



Fig. 3. Mass function. Solid line: solution of the renormalised 
gap equation using using Eqs. (^|), (|lo|); data: lattice simula- 
tions |lq], m = 0.036/a ~ 60MeV; dashed-line: gap equation 
solution using Eqs. (|), (|l|) (In the DSE studies, the renor- 
malisation point C = 19 GeV and the current-quark mass is 
0.6m* GcV [Eq. ©].) 



2.4 Gap equation and QCD's gauge sector 

Contemporary Landau gauge DSE studies of QCD's gauge 
sector have been used to infer a form of the effective inter- 
action pj| , which is well represented by (x = <3 2 Mqcd) 



— Q{x) := a(x) = 



a + a±x 



1 + a,2X + a^x 2 + a^x 3 ln(e + x) ' 

(14) 



an 



a i 



«2 



f'3 



«4 



1.47 0.881 0.314 0.00986 0.00168 



(15) 



Aspects of these DSE studies are in 
with quenched lattice-QCD results 



ualitative agreement 
M; e.g., the feature 



that the dressed-gluon propagator is finite at Q = 0. This 
effective interaction yields the results in the second row of 
Table ^: apparently, f„ is too small by a factor of ~ 2 and 
the scale of DCSB too small by a factor of ~ 4 = 1.6 3 pij . 
This outcome could have been anticipated from Ref. [[7L 

Equations (^) , (||) are correct on the ultraviolet domain 
in which perturbation theory is applicable. However, mod- 
elling is involved in extrapolating onto the infrared do- 
main. One may require that multiplicative renormalisabil- 
ity of the dressed fermion propagator be preserved in this 
modelling |22j . The Ansatz of Ref. |2^] incorporates this 
constraint, however, it implicitly specifies a kernel in the 
gap equation that violates a more important constraint, 
namely, Lorentz covariance [ p4[ . A minimal Ansatz that 
satisfies both constraints and Eq. (0) finds expression in 
the kernel of the renormalised gap equation as 



47ra(Q 2 ) 



1 



free 



Z 2 {Q 2 ) D ^ 



<Q)y7^(<?)y7,- (16) 



The gap equation's solution now gives the results in the 
third row of Table |[ [We simplified the numerical study by 
using a well-tried angle approximation: a((k—q) 2 )/Z 2 ((k— 
q) 2 ) k, a{{k — q) 2 )/Z 2 (ma,x(k 2 ,q 2 )).] This procedure is 
evidently a candidate for providing a bridge between DSE 
studies of the gauge sector and what is known, from the 
calculation of observables, to be the infrared strength re- 
quired in the effective interaction. The increased interac- 
tion tension arises because the enhancement of 1/Z(p 2 ), 
evident in Fig. ||, persists with this kernel and amplifies 
the effective coupling: the interaction tension calculated 
with Eq. ( |i"4| ) alone is o A = 12 GeV/fm but using Eq. 
([lq) we estimate that the l/Z 2 (Q 2 ) factor raises this to 
cX~25GeV/fm. 



3 Confinement 

This is a statement about the properties of coloured Sch- 
winger functions; i.e., gauge-dependent quantities. Hence 
it is plausible that the expression of confinement will de- 
pend on the gauge fixing procedure one employs. Linear 
covariant gauges are most widely employed in DSE stud- 
ies. While the problem of Gribov copies is unresolved, it af- 
fects only the far infrared behaviour of coloured Schwinger 
functions and hence is hoped to have little impact on phys- 
ical observables. This conjecture is supported by the close 
correspondence between Landau and Laplacian gauge two- 
point quark and gluon Schwinger functions obtained in 
quenched lattice-QCD simulations (^5). (Laplacian gauge 
fixing is a nonlocal prescription that is identical to Landau 
gauge at ultraviolet momenta. It is free of Gribov copies.) 

It has long been known |26| that for confinement it 
is sufficient that no coloured Schwinger function possess 
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a spectral representation. This is equivalent to the state- 
ment that all coloured Schwinger functions violate reflec- 
tion positivity. That property can manifest itself almost 
as simply as, say, dressed-gluons being described by a 
propagator with a large width. The connection between 
confinement and the violation of reflection positivity has 
cleanly been illustrated in QED3 ^7j, and it is a result 
of DSE studies that a sufficiently large interaction tension 
will always produce a dressed-quark propagator that vi- 
olates reflection positivity. (NB. Gap equation solutions 
that violate reflection positivity also exhibit DCSB. The 
converse is not necessarily true.) Thus it is a widespread 
conjecture that confinement in QCD is realised this way. 

A Schwinger function that violates reflection positiv- 
ity is easily identified: it is sufficient that the function 
exhibit a maximum or inflection at some spacelike four 
momentum-squared. In lattice-QCD simulations, this fea- 
ture is evident in the two-point gluon Schwinger function 
po| and incipient in the scalar functions characterising 
the dressed-quark two-point function Jl6[ |. It is plain in 
our results for the dressed-quark propagator, Figs. 0, [| 
and apparent in the DSE result for the gluon two-point 
function (^8). The DSE solution for the ghost two-point 
function also violates reflection positivity but in this case 
it is manifest in the strength of its singularity at k 2 = 0. 

A consistent picture may be emerging. It is comple- 
mented by the fact that the colour-singlet S'-matrix el- 
ements which describe physical processes do not exhibit 
unphysical quark and gluon production thresholds when 
calculated using coloured Schwinger functions that vio- 
late reflection positivity. The ideas we have described also 
ensure that coloured composites such as diquarks do not 
appear in the spectrum P,pT[. 



4 Epilogue 

The momentum-dependent dressing of quark and gluon 
propagators is a fact. It is certainly the keystone of DCSB 
and quite likely plays a central role in confinement. Its 
incorporation into models of low-energy phenomena fa- 
cilitates the development of a qualitatively reliable un- 
derstanding and makes possible an excellent description 
of experiment. That has provided guidance for modelling 
the long-range behaviour of the interaction between light- 
quarks. A remaining challenge is to calculate it and pro- 
gress is being made. 
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